Monte-Carlo simulation of localization dynamics of excitons in ZnO and CdZnO quantum well 

structures 
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Localization dynamics of excitons was studied for ZnO/MgZnO and CdZnO/MgZnO quantum wells (QW). 
The experimental photoluminescence (PL) and absorption data were compared with the results of Monte Carlo 
simulation in which the excitonic hopping was modeled. The temperature-dependent PL linewidth and Stokes 
shift were found to be in a qualitatively reasonable agreement with the hopping model, with accounting for 
an additional inhomogeneous broadening for the case of linewidth. The density of localized states used in the 
simulation for the CdZnO QW was consistent with the absorption spectrum taken at 5 K. 



Recently, there have been extensive studies on optical prop- 
erties of semiconductor quantum structures due to their size- 
dependent optical responses |lfl. In particular, wide-gap II- VI 
oxides have attracted a great deal of attention owing to the fact 
that their excitons have large binding energies |2LuD- The re- 
cent success of p-type doping [2, 3] in ZnO has opened a door 
to realize practical application of blue and ultraviolet (UV) 
light-emitting diodes, which can be an alternative to those 
based on Ill-nitrides. Introduction of cadmium and/or mag- 
nesium into ZnO plays a key role in strain ]4j] and band engi- 
neering of these oxides 1511a]. In previous work, we have re- 
ported an "anomalous" S-shaped temperature dependence of 
the photoluminescence (PL) band peak that has been observed 
in ZnO/MgZnO and CdZnO/MgZnO multiple quantum wells 

(mqw) am. 

Here, in addition to the S-shaped dependence, we report the 
observation of a W-shaped temperature dependence of the PL 
linewidth 1 10] . For quantitative discussion, we model these 
anomalous temperature behaviors of the emission band and 
compare the results with the experimental data. We use a 
Monte Carlo simulation which is based on the phenomeno- 
logical treatment of the excitonic kinetics 1 1011. The compari- 
son between simulation results and experimental data enables 
us to deduce the energy scale of distribution of the localized 
excitonic states. 

The samples investigated here were ZnO and CdZnO QWs, 
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whose growth and characteristics have been described else- 
where. Both of these QWs were grown with laser molecular- 
beam epitaxy on ScAlMg04 substrates, while the CdZnO QW 
includes a ZnO buffer layer. The material of our barrier layers 
is 5-nm-thick Mg x Zni_ x O. The PL was studied in the tem- 
perature range from 5 to 300 K. The continuous-wave He-Cd 
laser (emitting wavelength being 325 nm) was used full . 



The 



luminescence spectra taken for 



the 



ZnO/Mgo.27Zno.730 QW (the well width being 1.8 nm) 



(the 



exhibit a double band (8J, 11 211 at intermediate temperatures 
of 95 to 175 K. Figure 1 shows a typical near-band edge 
PL spectrum measured at 120 K. These two emission bands 
correspond to radiative recombination of localized excitons 
(LXi and LX2). The temperature dependence of the peak 
(see also Fig. 2 of Ref. 0) exhibits a well-known S-shaped 
variation 1 8], it slightly red-shifts in the range from 5 to 50 K, 
then blue-shifts in the range up to «100 K, splitts into a 
doubled peak at temperatures of 100 to w200 K, and then 
red-shifts again at higher temperatures. The initial red-shift 
reflects the ability of localized excitons to reach lower-energy 
sites via thermally activated hopping. The blue shift and 
the splitting of peaks that are observed with further increase 
in temperature may be attributed to the population of the 
localized exciton states with increasingly higher energy, 
whereas the subsequent red-shift is due to the typical band 
gap shrinkage with temperature. 

The observed features of hopping of localized excitons in 
our QWs are supported by the line width and Stokes shift 
measurements. The points in Fig. 2(a) and 2(b) show temper- 
ature variations of the full width at half maximum (FWHM) 
and Stokes shift of the PL band. Here, we notice a W- 
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FIG. 1: The excitonic absorption and PL spectra in a 
ZnO/Mgo.27Zno.730 MQW (well width being 1.8 nm) taken at 
120 K. 



shaped (decrease-increase-decrease-increase) temperature de- 
pendence 1 10] of the line width with a characteristic plateau at 
temperatures of 100 to 175 K. Its abrupt increase with temper- 
ature (5-100 K) was related to thermalization of the excitons 
over localized states, as opposed to the situation at lower tem- 
peratures where the excitons do not attain their equilibrium 
distribution. This is a typical manifestation of hopping of lo- 
calized excitons. On the other hand, the temperature variation 
of Stokes shift is more complicated due to the peak splitting 
at the intermediate temperatures. We show that these exper- 
imental results qualitatively agree with Monte-Carlo simula- 
tions explained later. 

The observed temperature dependences were simulated by 
a two-dimensional (2D) Monte Carlo algorithm. This algo- 
rithm adopts the Miller-Abrahams rate for phonon-assisted 
exciton tunneling between initial and final states (/ and j) with 
the energies of q and Cj , respectively ifToll : 



vq exp 



2r, ; 



2kT 



(1) 



Here is the distance between the localization sites, a is 
the decay length of the exciton wave function, and vq is the 
attempt-escape frequency. Hopping was simulated over a ran- 
domly generated set of localized states with the sheet density 
of N. Density of states (DOS) of the localization energies was 
assumed to be in accordance with a Gaussian distribution, 



g{e) = (N 2 /2na 2 ) 2 exp (-(e - E ) 2 /2a 2 ) 



(2) 



with the peak positioned at the mean excitonic energy Eq and 
the dispersion parameter (the energy scale of the band poten- 
tial profile fluctuation) a. The reference energy was e = Eq, 
which is the center of the distribution. All the energies in our 
simulation were below this value; i.e., we employed a "half" 
distribution II 131. In the case of the Gaussian DOS, the en- 
ergies were also assumed to be uncorrelated. For each gen- 
erated exciton, the hopping process terminates by recombina- 
tion with the probability Tq 1 and the energy of the localized 
state, where the recombination has taken place, is scored to 
the emission spectrum. 
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FIG. 2: (a) Evolution of the FWHM of the PL band with tempera- 
ture in a ZnO MQW. Open squares show experimental data; dashed, 
dotted, and solid lines depict results of the Monte Carlo simulation 
of the exciton hopping for different scales of random potential fluc- 
tuations a (indicated) with the inhomogeneous broadening Y taken 
into account; the dash-dotted line represents results for a=20 meV 
and r=0. (b) Simulated and experimental temperature dependences 
of the Stokes shifts of the PL positions with respect to the excitonic 
transition energies. 
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FIG. 3: Same as Fig. 2 except the sample is now 
Cdo.04Zno.96O/Mgo.12Zno.8sO QW (well width being 1.8 nm). 



The temperature variation of the data is basically deter- 
mined by the spatial (the product A^a 2 ) and temporal (the 
product tqVq) parameters of the hopping processes. For ex- 
ample, the kink in the temperature dependence of the FWHM 
is related to the energy of the potential fluctuations, a, i.e. 
a = 2fce7kink (Ref. fiol fl4l) . A slightly decreasing behav- 
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ior just above the saturation temperature (Tkink) indicates the 
formation of a thermalized exciton distribution. 

It has been difficult to perfectly reproduce both of our 
FWHM and energy shifting bahaviors with the identical set of 
simulation parameters. Nevertheless, if we adopt the values of 
Na 2 — 1, t i> = 2xl0 3 , and a— 20 meV, the calculated result 
approximately reproduces the 22-meV amplitude of the vari- 
ation of FWHM curve as shown by a solid curve in Fig. 2(a). 
The dashed and dotted lines in Fig. 2(a) demonstrate the sen- 
sitivity of the simulated results with respect to the hopping 
energy scale a. Increased a matches better with the experi- 
ment, but, as explained later, we will meet serious difficulty in 
accommodating the experimental Stokes shifts. The quanti- 
tative agreement with the experimental data requires an addi- 
tional inhomogeneous broadening (T = 12 meV) to be intro- 
duced 1 10]. In other words, our Monte Carlo simulation anal- 
ysis reveals that two kinds of inhomogeneities (a and T) are 
inherent to our QW. The dash-dotted line shows the simulated 
dependence without inhomogeneous broadening (V). The fur- 
ther increase of the linewidth above 175 K is attributed to the 
participation of the longitudinal optical phonons in the radia- 
tive transition and to the influence of delocalized excitons that 
are not accounted for in the model adopted here. 

The temperature behavior of the Stokes shift deduced from 
the peak position of the simulated spectra is depicted in 
Fig. 2(b). The following feature was reproduced in the sim- 
ulated results; the crossover temperature for the energy min- 
imum in the S-shaped dependence (50-60 K). On the other 
hand, the initial red-shifting behavior at low temperatures is 
somewhat exaggerated. We can accomodate this disagree- 
ment if the value of tqVq is reduced. But, in such a case, the 
agreement for the FWHM data becomes worse. Moreover, 
the simulation results do not reproduce the peak splitting phe- 
nomenon observed at the intermediate temperatures, the curve 
of which approximately runs in the half way of the dual peaks 
(open circles and squares). The more refinement of the model 
should be necessary to describe this splitting. Nevertheless, 
it is thought that the simulated Stokes shift energies corre- 
spond to the average energy of the splitted peaks observed 
in the experiments. At temperatures higher than 230 K, the 
experimental data of the Stokes shift are not in good agree- 
ment with those of simulation. This is probably because the 
spectral weight moves to the direction of the emission peak 
of longitudinal-optical phonon replicated recombination at the 
elevated temperatures as has been pointed out for ZnO bulk 
crystals and thin films 

For comparison, the simulated results for the 
Cdo.04Zno.96O/Mgo.12Zno.88O QW (the well width be- 
ing 1.8 nm) were shown in Fig. 3(a) and 3(b). The solid 
line in Fig. 3(a) represents the best fit obtained for the 
following values of the parameters: Na 2 =l, to^o = 7xl0 3 , 
r = 29 meV, and <r=20 meV. Unfortunately, the experimental 
Stokes shift datum was obtained solely at 5 K. This is in a 
good agreement with the result of simulation as shown with 
the open circle in Fig. 3(b). 

Finally, Fig. 4 compares the 5 K absorption spectrum (a 
dashed line) with the actual density of localized states used in 
the simulation (a solid line), 



CdZnO/MgZnO QW , \„„„« 




3.0 3.1 3.2 3.3 3.4 

Photon energy (eV) 



FIG. 4: Optical absorption spectrum (dashed line) of a CdZnO MQW 
taken at 5 K. A solid line shows the density of band-tail states used 
in the simulation with the inhomogeneous broadening (r) accounted 
for. The most prominent peaks ("ZnO A, B") come from the ZnO 
buffer layer. 
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where the inhomogeneous broadening is taken into account by 
introduction of T (here 29 meV). These densities of band-tail 
states are centered at the exciton transition energy (3.285 eV) 
obtained by adding the simulated Stokes shift [Fig. 3(b)] with 
the experimental PL energy (see also Fig. 2 of Ref. 7). In ad- 
dition, the long-wavelength wing of D(E) is seen to be in a 
fairly good coincidence with the low-energy tail of the corre- 
sponding absorption spectrum. This is thought to be an addi- 
tional proof on the validity of the exciton hopping model. 

The value of T for the CdZnO QW (29 meV) is significantly 
larger than that for the ZnO QW (12 meV), which is not sur- 
prising because of the more severe inhomogeneity in the for- 
mer sample; both the well and barrier layers are comprised of 
alloyed materials. Accordingly, to^o for the CdZnO QW is 
also larger. It is likely that the effective lifetime (related to To) 
is longer because it is known that the non-radiative decay path 
becomes reduced for more deeply localized excitons 0. 

Here we try to evaluate precise Cd concentration at irradi- 
ated region of the sample from a variational calculation fl7Tl 
where the electron-phonon interaction is properly taken into 
account. In general, the precise determination is difficult for 
Cd doped samples because of the reevaporation during the 
growth or to the large inhomogeneity. If the concentration 
is set to 4.3%, the transition energy is calculated to be 3.288 
eV, which is very close to the exciton transition energy at 5 K. 

In conclusion, we have studied evolution of the PL maxi- 
mum and FWHM of ZnO and CdZnO QWs with temperature 
by the Monte Carlo simulation. By making efforts to repro- 
duce both of our experimental FWHM and energy shift data 
with the same set of parameters consistently, characteristic en- 
ergy scale of the distribution of the localized states for the 
ZnO QW was evaluated to be a w20 meV. The simulation 
revealed two kind of inhomogeneities (a and T) are inherent 
to our ZnO-related QWs. The values of T and to^o for the 
CdZnO QW are significantly larger than those for the ZnO 
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QW. The density of localized states employed in the simu- 
lation turned out to be in a good agreement with the long- 
wavelength region of the absorption spectrum. 
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